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Technical Field of the Invention 

The present invention relates to a channel 
estimator and, more particularly, to a channel estimator 
that estimates the response of a channel through which 
10 signals are received. 

Rackaround of the invention 

For communication systems utilizing a 
periodically transmitted training sequence, least-squares 
15 algorithms or correlation algorithms have been the most 
widely used alternatives for the estimation of the 
channel. The estimation of the channel may then be used 
for various purposes such as initializing the values of 
the taps of an equalizer at a receiver. 
20 Both algorithms use a stored copy of the known 

transmitted training sequence at the receiver. The 
properties and length of the training sequence generally 
depend on the particular standard specifications of the 
communication system. In addition, a training sequence 



10 



is also used in a communication system for 

synchronization to indicate the beginning and/or end of a 

transmitted data frame or packet. 

in general, a communication system 10 that 

utilizes an equalizer adjusted by use of the training 
signal is shown in Figure 1 and includes an equalizer 12 
and a channel estimator 14 of a receiver. The base-band 
transmitted signal s(t) having a data rate of 1/T symbols 
per second is represented by the following equation: 



where I k e A -{«,.•••.«„> <= C is the transmitted training 
sequence, which is a discrete M-ary sequence taking 
15 values on the generally complex M-ary alphabet A , which 
also constitutes a two dimensional modulation 
constellation. The function q(t) as shown in a block 16 
represents the pulse shaping performed by a transmitter 
filter of finite support [-T q /2, T g /2] . 
20 The overall complex pulse shape is denoted by 

p(t) and is given by the following equation: 

p(t) = q(f)*q*(-t) {2) 



-2 
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where q\-t) shown in a block 18 is the receiver matched 
filter finite impulse response. Although it is not 
required, it is assumed, for the sake of simplifying the 
5 notation, that the span, T q , of the transmitter and 
receiver filters is an integer multiple of the symbol 
period T; that is, T q = N q T, N q &Z* . 

The physical channel between the transmitter 
and the receiver is denoted by c(f,T) and is shown by a 
10 block 20, and throughout the discussion below the 

concatenation of p(t) and the physical channel c(t,T) will 
be denoted by h(t) and is defined by the following 
equation: 



15 



h(t) = q{t) * c(r,T) *q (-0 = pit) * c(f ,r) < 3 ) 



where h(t) is generally referred to as the channel 
response. The physical channel c(t,T ) between the 
transmitter and the receiver is generally described by 
20 the following impulse response: 



c(*,t)= %c k <r)S(t-r k ) (4) 
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which describes a time -varying channel, where {c*(t)}cC 1 , 
-K<k<L, and t y re R . However, it is assumed that the 
time-variation of the channel is slow enough that 
c(r,r) = c(0 can be assumed to be a fixed (static) inter- 
5 symbol interference (ISI) channel throughout the training 
period; that is, it is assumed that c k (t) = c k , which in 
. turn implies that c(t) is defined by the following 
equation: 

10 c(0= £c*S(f-T*) < 5 

k=-K 

for 0<t<NT , where N is the number of training symbols. 

The summation indices K and L, where K is the 
number of anti-causal multi-path delays and L is the 
number of causal multi-path delays, are dependent on the 
15 channel and are more fully explained below in the 

discussion of correlation based channel estimation and 
subsequently. The multi-path delays r k are not assumed 
to be at integer multiples of the sampling period T. 
Indeed, an accurate and robust way of recovering the 
20 pulse shape back into the concatenated channel estimate 
when the multi-path delays are not exactly at the 
sampling instants is disclosed herein. 
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Equations (3) and (5) can be combined according 
to the following equation: 



L 

h(t) = p(t)*c(t)= YjC k p(t-T k ) . (6) 



5 



Because both p(t) and c(t) are complex functions, the 
overall channel impulse response h(t) is also complex. 
Therefore, p(t) , c(t), and h(t) can be written as 
functions of their real and imaginary parts as given by 
10 the following equations: 

p(t) = Pl (t) + jp Q {t) (7) 

c(0 = c/(0 + ^ G (0 (8) 

15 

Substituting equations (7) -(9) into equation (6) produces 
the following equation: 

h (t) = p{t) * c(t) - (/>/ (0 + JPq (')) * (c/ (0 + jc Q (0) 
20 =(MQ*C/(0- Pg(0*f£(0) + j (Pi(0*Cfl(0 + PflCO*g/(0) {10) 

' * *e w 
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The output matched filter 18 is denoted y(t) and is given 
by the following equation: 

y (t) = ^Lh^(.t - mym +v(o ( 11} 

5 

where v(J) designates the complex noise process after the 
matched filter. As in the case of equations (7)- (9), 
y(t) can be written in terms of its real and imaginary 
parts as y<0« y,(0 + /y G (0 • Sampling the output of the 
10 matched filter 18 at the symbol rate using a sampler 22 
produces the discrete-time representation of the overall 
communication system according to the following equation: 

*„]. y(0Ur ="Zhh[n-k]+v[n] . (12) 

it 

15 

The transmitted symbols are composed of frames 
(or packets) of length n\ where the first N symbols are 
the training symbols a k . Within a frame of length »\ the 
symbols may be denoted by the following equation: 

20 



U. = 



a k , forO<k<N-l (13) 
d k , for N <n< N'-l 



-6- 
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where the distinction of the first N symbols is made to 
indicate that they are the known training symbols, and 
where it is possible that the training symbols a k belong 
to a subset of the M-ary constellation alphabet A; that 

5 is KSACA^O, «„)}• m fact, for the 8-VSB digital 

television system, the signal alphabet is A s {±1,±3,±5,±7} , 
while the training sequence can only take binary values 

within the subset A s {-5, + 5} . 

Figure 2 shows a least- squares channel 
10 estimator 30 that has been used for the channel estimator 
14 of Figure 1 . In order to fully estimate the channel 
defined by equation (6), the least-squares based channel 
estimation algorithm assumes that the starting point and 
the ending point of the channel taps are either known or 
15 can be bounded. This assumption plays a critical role in 
the overall quality and robustness of the least-squares 
estimation procedure. It can be equivalents assumed 
that the impulse responses of the discrete- time 
equivalent composite channel h[n] can be written as a 
20 finite dimensional vector according to the following 
equation : 

h=[h[-N a ],h[-N a+ l] h[-l],h[0],h[l] h[N c -l],h[N c ]] T 

(14) 



10 



where N a denotes the number of anti-causal taps of the 
channel, N c denotes the number of causal taps of the 
channel, and N a+ N c +1 is the total memory of the channel 

Based on equation (12) and assuming that 
N>N a +N c +\ where N again is the number of training 
symbols, the output of the matched filter 18 
corresponding to the known training symbols can be 
written as follows: 

y[N c ] = h[-N a ]a Nc+Nc +- + h[0]a Nc + - + h[N c }a Q +0[iV c ] 
y[N c +l] = h[-N a ]a Nc+Nc+1 + - + h[0]a Nc+l+ - + h[N c ] ai + v[N c +1] 



y[N 



15 



20 



which can be compactly written according to the following 
vector equation: 

y = Ah + 0 < 15) 



-8- 
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where y comprises values extracted by a sample extractor 
32 from the received signal and is defined according to 
the following equation: 



y = [y[N c ], y[N c + 1],. • • . y\. N ~ 1 - N a H 7 



(16) 



where A is based on the training signal and is defined 
according to the following equation: 



10 



A = 



a N c +N a a N c +N a -l 
a N c +N a +l a N c +N a 



a N-l 



a N-2 



a 0 

a-N-\-N a -N c ] iN - Na -N c )x(N a +N c +l) 



(17) 



and where v is defined according to the following 
equation: 



15 



v =[v[N c ]MN c +1] v[N-l-N a ]f 



As long as the matrix A is a tall matrix and of full 
column rank such that N>2(N a+ N c ) + l and rank{A}=N a+ N c+ l , 
then the least-squares solution which minimizes the 
20 objective function given by the following equation 



-9- 



^(h) = «y-Ah|| 2 



(18) 



exists and is unique and is given by the following 
equation: 



5 



h^=(A w A)- 1 A w y (19) 



where the superscript H denotes Hermitian transpose, and 
where is the least- squares estimation of the channel 

10 between the transmitter and the receiver. Thus, in . 
accordance with equation (19), a matrix multiplier 34 
multiplies the matrix y by a matrix A + , which is (A" A)" 1 A" 
and which is stored in a memory 36. Thus, the matrix 
multiplier produces the channel estimate h u according to 

15 equation (19) . 

One of the problems associated with the least- 
squares based channel estimation is the uncertainty 
associated with the actual channel spread, i.e., the 
spread between the outermost mul ti -paths . That is, this 
20 uncertainty arises because the channel memory and the 

first and last channel taps are not actually known by the 
receiver. Therefore, in order to set up equations (16)- 
(19), some sort of preprocessing is needed, such as 



10 
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correlating the stored training symbols with the received 
signal, or some pre-assumed starting and ending points 
for the channel must be used such that the total channel 
spread is not too long that would otherwise result in the 
5 A matrix being a wide matrix rather than a tall one. 

Even if the starting and ending points of the 
actual channel spread are known by the receiver, least- 
squares based channel estimators still have problems that 
may be insurmountable. For example, least- squares 
10 channel estimators cannot properly estimate the channel 
from the training sequence if the channel spread is too 
long for the training sequence to effectively support the 
channel . 

That is, if N < 2(N a +N c )+l, such that the number 
15 of training symbols is insufficient to support, or 

create, an over-determined channel convolution matrix A, 
then there are two cases to consider. 

First, if the full A matrix as given by 
equation (17) is used, the least-squares estimate of 
20 equation (19) no longer exists. That is, the matrix A H A 

is singular. 

Second, if the receiver does not attempt to 
construct the full A matrix and instead constructs an 
alternative thinner matrix A with N > 2(N a + N C ) + 1 such 



that N a <N a and/or N C <N C , then the inverse of A H A will 
exist but the new least-squares solution given by the 
following equation: 

5 h w =(A // A)" 1 A i/ y 

will not be a very desirable estimate. That is, the 
resulting channel estimate no longer resembles the actual 
channel, especially if the actual channel spread is much 
10 longer that the training sequence can support. 

An example is illustrated in Figures 3 (a) -3(e). 
A channel with delays T t = [-105.32. -87.643, -45.443, 0, 
12.22, 39.9, 203.659, 321.82, 435.12] with a complex 
raised cosine pulse shape of 41 samples (a channel spread 
15 of 580 symbols including the tails of the pulse shape) 
was simulated. The real part of the channel impulse 
response seen at the output of the matched filter is 
shown in Figure 3(a) . As shown in Figure 3(a), each peak 
of the channel impulse response has a raised cosine pulse 
20 shape. The training sequence is the same as in the 8 VSB 
digital television standard and is of length N = 704 
symbols. With 704 training symbols, a channel with a 
maximum delay spread of only N a + N c + 1 = 352 taps, 
including the tails of the pulse shape, can be supported. 

- 12- 
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Thus, the actual channel has a much longer delay spread 
than the standard least-squares based method can handle. 
The least-squares channel estimate corresponding to cases 
where the receiver constructs a smaller convolution 
5 matrix A with -AT, =-125, # c =60, -# fl =-125, N c =120 , 

_N a =-125, # c =206, and -N a =-125 , tf c =342 are shown in 
Figure 3(b), Figure 3(c), Figure 3(d), and Figure 3(e), 
respectively. The matrix A is non-singular for all cases 
except Figure 3(e) . As can be seen from Figure 3(b) to 
10 Figure 3(d), as the assumed channel spread gets longer, 

the estimated channel response increasingly deviates from 
the original channel response. In Figure 3(e), where A 
is undetermined, the estimated channel no longer 
resembles the original channel response. However, it may 
15 be misleading to assume a very short channel spread in 
order to obtain a cleaner channel estimate within the 
assumed channel span. In that case, it is very likely 
that very significant early pre-cursor multi-paths or 
late post-cursor multi-paths may be missing. 
20 as indicated above, it is also known to apply 

auto-correlation in estimating the channel. The sampled 
received sequence at the output of the matched filter is 
denoted by y[n] and is given in equation (12), where the 
pre-cursor symbols are at positions «>0 and post-cursor 
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symbols are at positions n<0. The known training 

sequence is denoted by s [n] = a n f or n = 0 N-1. 

The auto-correlation functions r s [m]and r sy [m] are defined 
in accordance with the following equations: 



N-l 

r s [m]= + m] 

N-1 
k=0 



(20) 



(21) 



10 



where {s[fc]e A} . 

It is assumed that, in order to be able to use 
correlation based channel estimation, the training 
sequence must belong to a certain class of sequences and 
thereby possess certain desirable properties. An 
exemplary class of sequences is maximal length pseudo- 
15 noise (PN) sequences. A PN-sequence of length n is 

denoted herein as PN n . In general, the periodic auto- 
correlation of a binary valued (f+A, -A}) PN n sequence is 
given by the following equation: 



20 



r PNn [m]= ■ 



A 2 n, ifm = 0,±n,±2n,. 
- A 1 , otherwise 



(22) 



-14- 
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Figure 4 shows a correlation channel estimator 
40 that can be used for the channel estimator 14 to 
estimate the channel in order to initialize the equalizer 
12. The correlation channel estimator 40 includes a 
5 cross-correlator 42 that cross-correlates the received 
signal with the training sequence stored in a memory 44 
and a level thresholder 46 that sets to zero any 
correlation coefficients that are below a predetermined 
threshold. Accordingly, the correlation channel 
10 estimator 40 produces a channel estimate h . 

More specifically, based on the auto- 
correlation property of equation (22) , the channel can be 
estimated by cross-correlating s [n] (the training 
sequence known and stored in the receiver) with the 
15 received sequence y[n]. The initial pre-cursor impulse 
response estimate h a [n] is determined from the cross- 
correlation of the training sequence s [n] stored in the 
receiver and the actual received symbols at lags from 0 
to -N+l with respect to the start of the received 
20 training sequence. That is, h a [n\ is given by the 
following equation: 



N-l 



h a [n] = -7- 5>* [fc]** + n] far lags n = 0,-1,... -N + 1 



(23 
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Similarly, the initial post-cursor response ■ estimate h c [n\ 
is determined from the cross-correlation of the training 
sequence stored in the receiver and the actual received 
5 symbols at lags from 1 to N-l with respect to the start 
of the received training sequence. That is, h c [n] is 
given by the following equation: 



,7 fn] = _!_ V S [k]y* [k + n] for lags n = 1,2..., N -1 
° U J r s [0] ^ 



(24) 



10 

where r 8 [0] in equations (23) and (24) is defined 
according to the following equation: 



N-l 2 



r s [0]=XlK| 



15 



Equations (23) and (24) imply that the correlations 
should be extended on either side by the length of the 
training sequence in order to obtain the longest possible 
channel estimate within the limits of the transmission 
20 standard, with the assumption that NZK and N>L. 

With respect to the channel estimation 
procedure given by equations (23) and (24), the overall 
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initial estimated channel transfer function H( Z ) is given 
by the Z-transform of the concatenated channel estimates 
cording to the following equation: 



ac 



= 2M*]z"* (25) 



where the concatenated estimated channel impulse fc[n] i- 
given by the following equation: 

10 

r r n _KM, M-iV+l<n<0 (26) 

Equation (25) can be rewritten in accordance with the 
following equation: 

15 

r 5 [0] 

where A is a (N + N a + N c ) by (N a + N c + D Toeplitz 
matrix with a first column [a 0 , a x , . • • . a N-i' °' ' ' " 
20 0] T and a first row [a 0 , 0; . . . . 0] T , and where y = 
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[y[-N a ] y[N + N c -1]] T . m order to get rid of 

the side lobes of the aperiodic autocorrelation, the 
autocorrelation of the training sequence is determined 
according to the following equation: 

i_A*A (28) 



10 



R„„ = — rrrA A 



aa r s [0] 



and then R aa can be inverted and applied to H u (z) 
according to the following equation: 



>" l w (*\ (29) 



H c (z) = R^H„(z) 



However, the channel estimate H c {z) obtained by equation 
(29) has the contributions due to unknown symbols prior 
15 to and after the training sequence, as well as the 

additive channel noise; only the side lobes due to the 
aperiodic correlation is removed. 

As discussed below, the correlations of 
equations (23) and (24) will clearly yield peaks at 
20 keZ+}, for k = -K, , - -1, 0, and at fceZ*}, for k 
= 1, . . L. These delays D a k and D{ are the sampling 
instants closest to the locations of the actual physical 
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channel multi-path delays X k for k = -K, . . ., -1. 0, 1, 
. . , L within a symbol interval . 

As shown in Figure 4, there is a level 
thresholding algorithm following right after the cross- 
correlations such that the estimated channel taps are set 
to zero if they are below a certain threshold. This 
thresholding algorithm operates in accordance with the 
following equation: 



10 seth c [n] = 



0, ifh c [n)<e (30) 
h c [n], otherwise 



, „ 1 n 1 .„ N-l, where the 

for n = -N+l, . - •/ - 1 ' °' x > ■ ' " 

subscript c denotes that the corresponding function has 
been cleansed of side lobes. As a result, the estimated 
15 transfer function H c {z) after thresholding can, in 

general, be written according to the following equation: 



20 



where for l<i<L, and fc-fc/Oo *>r l<k<K. 

The notation of equation (31) is adopted because, after 
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thresholding, there are usually very few dominant taps. 
Accordingly, the channel after correlation and 
thresholding is sparse. Moreover, the peaks in the 
channel estimate no longer have the pulse shapes of the 
5 actual channel response. 

in equation (31) , the coefficients P k ^ 
correspond to the post-cursor (causal) part of the 
correlation, and the coefficients a, correspond to the 
pre-cursor (anti-causal) part of the correlation. The 
10 delays Df are the anti-causal delays, and the delays Df 
are the causal delays. As a convention, it is assumed 
that 1 < D{ < Dl < • • • < D C L and simi larly that 
l<Df <D\ < -<D a K . Also, by the construction of the 
channel estimates in equations (23) and (24), D a K <N-l 

15 and D c L <N-l. 

At this point, the causal part H cc and the 
anti-causal part H ac of the transfer function of equation 
(31) can be defined in accordance with the following 
equations: 



20 



A.w-a.tf,*'* **** + » <32) 
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H cc ( Z ) = « 0 («i^ lC +a 2Z~ Dl +- + cc L z- D h (33) 



such that H c (z) = H ac (z) + H cc (z)- Again the first subscript 
denotes whether the part is anti-causal or causal, and 
the second subscript denotes that the part has been 
cleansed of the side lobes. The partition of equations 
(32) and (33) enable the development of the appropriate 
feed-forward filter initializations for a Decision 
Feedback Equalizer (DFE) . 

Unfortunately, if the PN training sequence is 
finite and the standard linear correlation is used, the 
auto-correlation values corresponding to the non-zero 
lags will not be constant and will not be as low as -A 2 . 
As a simple illustration, a sequence composed of six PN 51 i 
15 appended back to back may be considered so that the 
sequence is given by the following equation: 



10 



y = [/W 511 , PN 5U , PN 5U , PN m , PN 5ll ,PN sn ] . (34) 

20 Then, rxy [m], with x=[PN 511 ] T , will be given as shown in 
Figure 5. A low correlation value of -A 2 will be 
obtained for lags that are not multiples of L= 511 
corresponding to the intermediate PN 511 portions of the of 
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the long sequence y. However, as illustrated in Figure 
5, for the outermost lags, this constant and low 
correlation value of -A 2 will not be obtained. Instead, 
these correlation tails are noise like due to the 

5 finiteness of the training sequence. That is, the stored 
training sequence is correlated not only with the 
training sequence in the received signal but also with 
random data symbols and reserved symbols. 

The training sequence that is used for digital 

10 television is given by the following equation: 



S = [ss,PN 5U , PN 63 , ± PN & , PN 63 ] T 



where ss = [+5 -5 -5 +5] represents the four segment sync 
15 symbols at the beginning of the frame sync segment. 

There are also reserved frame bits and data bits right 
before and after the training sequence S . 

In summary, the correlations of the received 
signal with the stored training sequence will be noisy 
20 because the PN sequences are finite in length so that 

they will not achieve the desired low correlation values 
for non-zero lags, and because the span of the cross- 
correlation includes the known training sequence as well 
as random data symbols and reserved symbols . 
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Because of these reasons, the initial estimate 
of the channel that is acquired from the correlation of 
the matched filter output with the stored training 
sequence will not yield a clean and accurate channel 
5 estimate. Accordingly, even if the channel tap estimates 
R,&} are scaled properly, they will not be very close 
to the actual channel tap values {cnpCO)} due to the 
reasons given above as well as due to the fact that there , 
may also be tap values in the tails of the adjacent 
10 multi-paths. Hence, both the amplitudes and phases of 
the estimated tap values may be incorrect. 

Moreover, after the thresholding takes place, 
the most likely result is that the thresholded 
correlation will include only the peaks of the channel 
15 taps, which implies that the tails of the raised cosine 
pulse shape that are under the noisy correlation output 
may be eliminated entirely. This elimination of the 
tails may lead to critical performance loss in the 
initialization of the equalizer 12 and a slower 
20 convergence onto the correct tap values for the equalizer 
12 after initialization. 

The present invention mitigates the problems of 
the least-squares channel estimators and the correlation 
channel estimators providing a channel estimator that is 
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a blend of a least- squares channel estimator and a 
correlation channel estimator. 



Summary of the Invention 

5 In accordance with one aspect of the present 

invention, a method for determining a channel response 
comprises the following: determining a matrix T based on 
a correlation of a received signal and a known training 
sequence; extracting a vector y from a received signal; 

10 and, estimating the channel response from a least-squares 
solution based on the matrix T , the vector y, and a 
matrix A formed from the elements of the known training 
sequence . 

In accordance with another aspect of the 
15 present invention, a method for determining a channel 

response comprises the following: determining a matrix T 
based on a correlation of a received signal and a known 
training sequence; extracting a vector y from a received 
signal; determining an unknown coefficient vector Yls 
20 according to the following equation: 



Yls 



= (r H A H Ar)T 1 r w A // y 
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wherein the unknown coefficient vector Yls comprises a 
least-squares solution of the following equation: 

y = Ary+v 

5 

and wherein A comprises a matrix based on the training 
sequence; and, determining a channel response according 
to the following equation: 

!Q Kew = r YLS ■ 

in accordance with still another aspect of the 
present invention, a method for determining a channel 
response comprises the following: correlating a received 
signal with a training sequence; constructing a matrix T 
15 by approximating a desired pulse shape for each peak of 
the correlation; constructing a matrix A based on the 
training sequence; solving a least-squares solution of 
the following equation in order to derive an unknown 
coefficient vector y : 

20 

y = Ary 
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wherein y is a vector representing a received signal; 
and, determining the channel response from the unknown 
coefficient vector y and the matrix T . 

5 Brief Description of th e Drawings 

These and other features and advantages will 
become more apparent from a detailed consideration of the 
invention when taken in conjunction with the drawings in 
which: 

10 Figure 1 illustrates a generalized system whose 

equalizer is initialized by a channel estimator; 

Figure 2 illustrates a least-squares channel 
estimator that has been used for the channel estimator of 
Figure 1; 

15 Figures 3 (a) -3(e) illustrate an exemplary 

simulation useful in explaining problems associated with 
the least-squares channel estimator of Figure 2; 

Figure 4 illustrates a correlation channel 
estimator that has been used for the channel estimator of 

20 Figure 1; 

Figure 5 illustrates an exemplary simulation 
useful in explaining problems associated with the 
correlation channel estimator of Figure 4; and, 
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Figure 6 illustrates a channel estimator 
according to an embodiment of the present invention that 
is a blend of a least-squares channel estimator and a 
correlation channel estimator. 

5 

Detailed Description 

Figure 6 illustrates a blended channel 
estimator 50 having a sample extractor 52, a matrix 
multiplier 54, a training sequence memory 56 that stores 
10 a training sequence, a cross-correlator 58, a level 

thresholder 60, a combiner unit 62, a multiplier 64, and 
a tap weight initializer 66. The sample extractor 52 
extracts values y from the output of a receiver's matched 

filter, and the matrix multiplier 54 multiplies the 
15 extracted values y by a matrix (AT) 0 , where the matrix 

(AT) 0 is calculated by the combining unit 62 in the manner 
described below. The cross-correlator 58 cross- 
correlates the received signal with the stored training 
sequence stored in the memory 56 and the resulting cross- 
20 correlation is thresholded by the level thresholder 60 so 
that the combining unit 62 can construct the matrix (Alf 
from the output of the level thresholdewr 60 and the 
training sequence stored in the memory 56. 
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Thus, the blended estimation procedure 
illustrated by Figure 6 starts from the exact same steps 
of the correlation based estimation as outlined by 
equations (23) and (24), and then the thresholding step 
5 is performed on the channel estimate as indicated by 
equation (27) in order to clean up the noise in the 
channel estimate. These steps produce the cleaned 
channel estimate H c (z) shown in Equation (31) . In 
equation (31) , . the value 1 denotes the main signal and 
10 the other values in equation (2) denote multi-paths. 
Equivalently, the cleaned channel impulse response 
estimate can be represented as shown in the following 
equation : 

15 h c =dc 0 [P K , 0^0, fe.-.ft, M O^AftL-saw, ^fi, a L 

Dfr-tf^l zeros Df-1 zeros Df-1 zeros o£-Z>£_ r l zeros 

(35) 



where the zeros represent the spaces between the peaks. 
20 The structure and notation of equations (31) and (35) are 
important to the derivation of the algorithm below. 

As indicated above, the pulse shape p(t) is 
convolved by the physical channel c(t) according to 
equation (6) . As also indicated above, this pulse shape 
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is lost due to the correlation and thresholding process. 
Therefore, it is desirable to recover this pulse shape 
for each of the multi-paths. Accordingly, the shifted 
and scaled copies of the pulse shape p(t) (shifted by x k 
and scaled by c k ) may be approximated by a linear 
combination of three pulse shape functions shifted by 
half a symbol interval (T/2) . This linear combination is 
given by the following equation: 



10 c k p(nT-t k )- 



-K<k<-1 

(36) 



where fa k) , K < k < lJ = _, C C 1 . The last part of equation (36) 
contains an approximation of the pulse shape of the 
transmitter/receiver filter combination. By making this 
15 approximation, the tails of the complex pulse shape p(t) 
can be efficiently recovered even though these tails are 
buried under the noisy output of the correlation 
' processing and are lost when uniform thresholding is 
applied. 

20 to arrive at this approximation, three vectors 

p k are introduced for k=-l, 0, +1. each containing T 



spaced samples of the complex pulse shape p(t) shifted by 
kT/2 (for k = -1, 0, +1) according to the following three 



equations : 



P-i = 



Po = [p(-w,r)-,p(-np(o)ip(r),-,pKr)f 



q 2 



Pi = 



p(-V-|)-,p(- r -|) P | 



(37) 



(38) 



(39) 



By concatenating these vectors, a (2N q+ l) x 3 matrix P is 
10 defined according to the following equation: 



P = fc>-i,Po>Pi] 



(40) 



This vector represents shifted and sampled values of a 

15 raised cosine pulse. 

A matrix T whose columns are composed of the 

shifted vectors Pk , where the shifts represent the 

relative delays of the multi-paths. This matrix T is 

given by the following equation: 

20 
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15 



r = 



0 



(D K +D L, )jc3 



0. 



( D L +D A'-l )jc3 



0 



D K* 3 



0 



D c L x3 



0 



(D£ + D£)*3 
0 r ^ P 

(41) 



where T is of dimension (d% + D c l +2N q +l)by. 3(K + L + l) 
5 dependent on the channel spread, and Omxn denotes an m by 

n zero matrix such as a D%-D%_ x by 3 matrix of zeros. 
The received signal vector y and the training sequence 
vector A are defined by equations (16) and (17), 
respectively. Because of the assumption that q(t) spans 
N q symbol durations, q[n] has N q + 1 sample points and 

p[n] in turn has 2N g + 1 samples. Therefore, N a = D% + N q 

and N c = D C L + N g . 

An unknown vector Y is defined in accordance 
with the following equation: 



10 



(42) 



where y ( *> = ,yf } -W^ \ for -K<k<L. The unknown vector 
Y has coefficients {yl*} where n = -1, 0, 1, and where k = 
_ K/ _ . .,o, . . . , L and has a length of 3 (K+L.+1) - 
Similarly to equation (15), the vector y representing the 
5 output of the receiver's matched filter may be given by 
the following equation: 

y = Ary+v < 43) - 

10 where v is the observed noise vector. Using a least- 
squares solution, the unknown coefficient vector y can be 
estimated in accordance with the following equation: 

y iS =(r«A w Ar)r 1 r // A w y = (Ar) 0 y (44) 

15 

Once the least-squares estimated vector Yls i« obtained, 
the new channel estimate h new , where the pulse tails are 
recovered, can simply be obtained in accordance with the 
following equation: 

20- 

f ro (45) 
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The multiplier 64 multiplies Yls from the matrix 
multiplier 54 by T according to equation (45) in order to 
produce h new . Equations (44) and (45) may be solved using 
either a singular value decomposition (SVD) method, by a 
5 conjugate gradient (CG) method, or by any other suitable 
method. 

Once the new estimated channel vector is 
determined per equation (45), the reconstructed vector at. 
the output of the receiver matched filter can be defined, 
10 based on h new according to the following equation: 

y = A/; iew (46) 
y = Ary^ (47) 

y^Ar^A^Ar^r^A^y (48) 



15 



The variance of the noise samples at the output of the 
receiver's matched filter, cr| , can be estimated based on 
the following equation: 



20 *i = - — -||y - yf ( « ) 



which arises naturally. The length of the vector Y is 

N-D a K -D c L . 

The tap weight initializer 66 may then 
initialize the tap weights of the taps of the equalizer 

5 12 using the channel estimate h new and the noise variance 

in a manner that is known in the art. 

When the present invention is used to 
initialize the taps weights of the taps of an equalizer, 
the noise variance given by equation (49) and the new 

10 channel response given by equation (45) can be combined 

in any known way so as to produce the initial tap weights 
for the taps of the equalizer. Alternatively, instead of 
combining the noise variance given by equation (49) and 
the new channel response given by equation (45), a small 

15 non-zero value noise variance may be assumed and this 
assumed value may be combined with the new channel 
response given by equation (45) so as to produce the 
initial tap weights for the taps of the equalizer. 

Modifications of the present invention will 

20 occur to those practicing in the art of the present 
invention. For example, the present invention as 
described above relies on the use of a training signal. 
However, other signals, such as data signals, can be used 
to implement the present invention. Also, although the 
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invention is specifically shown for use in initializing 
the taps of an equalizer, the channel estimate produced 
in accordance with the present invention can be used for 
other purposes . 
5 Accordingly, the description of the present 

invention is to be construed as illustrative only and is 
for the purpose of teaching those skilled in the art the 
best mode of carrying out the invention. The details may 
be varied substantially without departing from the spirit 
10 of the invention, and the exclusive use of all 

modifications which are within the scope of the appended 
claims is reserved. 
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We claim: 



1 1 . A method for determining a channel 

2 response comprising: 

3 determining a matrix T based on a correlation 

4 of a received signal and a known training sequence; 

5 extracting a vector y from a received signal; 

6 and, 

7 estimating the channel response from a least- 

8 squares solution based on the matrix T, the vector Y, and 

9 a matrix A formed from the elements of the known training 
10 sequence. 

1 2 . The method of claim 1 wherein the channel 

2 response has a pulse shape, and wherein the matrix T is 

3 constructed so as to approximate the pulse shape of the 

4 channel response. 
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1 3 . The method of claim 2 wherein the matrix T 

2 comprises three sets of samples derived from the pulse 

3 shape, wherein the samples have a period T, wherein the 

4 first set comprises the pulse shape samples shifted by 

5 +T/2, wherein the second set comprises the pulse shape 

6 samples not shifted, and wherein the third set comprises 

7 the pulse shape samples shifted by -T/2. 

1 4. The method of claim 3 wherein the pulse 

2 shape from which the samples are derived comprises a 

3 complex raised cosine pulse shape . 

1 5. The method of claim 1 wherein the matrix T 

2 is constructed with a number of zeros dependent on 

3 distances between peaks of the correlation. 

1 6. The method of claim 5 wherein the channel 

2 response has a pulse shape, and wherein the matrix T is 

3 constructed so as to approximate the pulse shape of the 

4 channel response. 
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1 7 . The method of claim 6 wherein the matrix T 

2 comprises three sets of samples derived from the pulse 

3 shape, wherein the samples have a period T, wherein the 

4 first set comprises the pulse shape samples shifted by 

5 +T/2, wherein the second set comprises the pulse shape 

6 samples not shifted, and wherein the third set comprises 

7 the pulse shape samples shifted by -T/2. 



1 8 . The method of claim 7 wherein the pulse 

2 shape from which the samples are derived comprises a 

3 complex raised cosine pulse shape. 

1 9 . The method of claim 1 wherein the least- 

2 squares solution is solved using a singular value 

3 decomposition method. 

1 10. The method of claim 1 wherein the least- 

2 squares solution is solved using a conjugate gradient 

3 method . 

1 11. The method of claim 1 further comprising 

2 determining initial tap weights for an equalizer based on 

3 the channel response. 
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1 12. The method of claim 11 further comprising 

2 determining the initial tap weights for the equalizer 

2 

3 based on the channel response and a noise variance 0*^ 

4 given by the following equation: 
5 

7 

8 wherein D denotes distance between peaks of the 

9 correlation, wherein the superscript a denotes a distance 

10 of a corresponding anti-causal peak, wherein the 

11 superscript c denotes a distance of a corresponding 

12 causal peak, wherein K denotes the number of anti-causal 

13 peaks, wherein L denotes the number of causal peaks, 

14 wherein N denotes a number of training symbols in the 

15 training sequence, and wherein y is based on the vector 

16 y. 
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1 13 . A method for determining a channel 

2 response comprising: 

3 determining a matrix T based on a correlation 

4 of a received signal and a known training sequence; 

5 extracting a vector y from a received signal; 

6 determining an unknown coefficient vector Yls 

7 according to the following equation: 
8 

9 Yls = ( r " A^ArjfY" A H y 

10 

11 wherein the unknown coefficient vector Yls comprises a 

12 least-squares solution of the following equation: 
13 

14 y = Ary + v 

15 

16 and wherein A comprises a matrix based on the training 

17 sequence; and, 

18 determining a channel response according to the 

19 following equation: 
20 

21 Kew =TYLS • 
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1 14. The method of claim 13 wherein the channel 

2 response has a pulse shape, and wherein the matrix T is 

3 constructed so as to approximate the pulse shape of the 

4 channel response. 

1 15. The method of claim 14 wherein the matrix 

2 T comprises three sets of samples p-i, p 0 / and p +1 derived 

3 from the pulse shape, wherein the samples have a period 

4 T, wherein the first set p_i comprises the pulse shape 

5 samples shifted by +T/2, wherein the second set p 0 

6 comprises the pulse shape samples not shifted, and 

7 wherein the third set p +i comprises the pulse shape 

8 samples shifted by -T/2. 

1 16. The method of claim 15 wherein the pulse 

2 shape comprises a complex raised cosine pulse shape. 

1 17 . The method of claim 15 wherein the matrix 

2 T comprises an approximated pulse shape P formed by 

3 concatenating p_i, Po/ and P+i • 
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1 



2 



18. The method of claim 17 wherein the matrix 
r is given by the following equation: 



5 
6 
7 



r = 



0 



0 



0 



(D L +D £-l )jc3 



0 

P 



0 



0 



D c r x3 



< D K + °£-l>* 3 
P 



0 



0 



P 



wherein the zero matrices are dependent on distances D 
between peaks of the correlation, wherein the superscript 

8 a denotes a distance of a corresponding anti-causal peak, 

9 wherein the superscript c denotes a distance of a 

10 corresponding causal peak, wherein K denotes the number 

11 of ant i- causal peaks, and wherein L denotes the number of 

12 causal peaks . 



1 19 . The method of claim 18 wherein the pulse 

2 shape comprises a complex raised cosine pulse shape. 
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1 



20. 



The method of claim 13 wherein the matrix 



2 



A is given by the following equation: 



3 



A = 



a N c +N a a N c +N a -l 
a N c +N a +l a N c +N a 



a 0 



4 



a N-\ 



a N-2 



a N-l-N a -N c 



(N-N a -N c )x(N a +N c +Y) 



5 

6 wherein ai comprise training symbols, wherein N comprises 

7 the number of training symbols, wherein N a comprises the 

8 number of anti-causal peaks in the correlation, and 

9 wherein N c comprises the number of causal peaks in the 
10 correlation. 

1 21. The method of claim 20 wherein the channel 

2 response has a pulse shape, and wherein the matrix T is 

3 constructed so as to approximate the pulse shape of the 

4 channel response . 
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1 



22. The method of claim 21 wherein the matrix 



2 r comprises three sets of samples p_i, p 0 , and p+i derived 

3 from the pulse shape, wherein the samples have a period 

4 T, wherein the first set p_i comprises the pulse shape 

5 samples shifted by +T/2, wherein the second set p 0 

6 comprises the pulse shape samples not shifted, and 

7 wherein the third set p +1 comprises the pulse shape 

8 samples shifted by -T/2 . 

1 23 . The method of claim 22 wherein the pulse 

2 shape comprises a complex raised cosine pulse shape. 

1 24. The method of claim 22 wherein the matrix 

2 T comprises an approximated pulse shape P formed by 

3 concatenating p_i, Po/ and p+i . 

1 25. The method of claim 24 wherein the matrix 

2 r is given by the following equation: 



3 



P 



0 



P 



0 



(Dj*+Df);c3 



0 




4 



r = 



0 



(D L +D A:-i )Jt3 



0 



<d£+d£_i>*3 
P 



0 



Djx3 



0 




0 



<°l- d £-i>* 3 
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5 

6 wherein the zero matrices are dependent on distances D 

7 between peaks of the correlation, wherein the superscript 

8 a denotes a distance of a corresponding anti-causal peak, 

9 wherein the superscript c denotes a distance of a 

10 corresponding causal peak, wherein K denotes the number 

11 of anti-causal peaks, and wherein L denotes the number of 

12 causal peaks. 

1 26. The method of claim 25 wherein the pulse 

2 shape comprises a complex raised cosine pulse shape. 

1 27. The method of claim 13 further comprising 

2 solving the equation for the unknown coefficient vector 

3 using a singular value decomposition method. 

1 28. The method of claim 13 further comprising 

2 solving the equation for the unknown coefficient vector 

3 /L^ using a conjugate gradient method. 

1 29. The method of claim 13 further comprising 

2 determining initial tap weights for an equalizer based on 

3 the channel response.. 
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1 30. The method of claim 13 further comprising 

2 determining the initial tap weights for the equalizer 

3 based on the channel response h new and a noise variance 

4 G$ given by the following equation: 
5 

6 <r|= lly - yf 

7 

8 wherein D denotes distance between peaks of the 

9 correlation, wherein the superscript a denotes a distance 

10 of a corresponding anti-causal peak, wherein the 

11 superscript c denotes a distance of a corresponding 

12 causal peak, wherein K denotes the number of anti-causal 

13 peaks, wherein L denotes the number of causal peaks, 

14 wherein N denotes a number of training symbols in the 

15 training sequence, and wherein y is based on the vector 

16 y. 
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1 4 31. A method for determining a channel 

2 response comprising: 

3 correlating a received signal with a training 

4 sequence; 

5 constructing a matrix T by approximating a 

6 desired pulse shape for each peak of the correlation; 

7 constructing a matrix A based on the training 

8 sequence; 

9 solving a least-squares solution of the 

10 following equation in order to derive an unknown 

11 coefficient vector y : 
12 

13 y = Ar> 
14 

15 wherein y is a vector representing a received signal; 

16 and, 

17 determining the channel response from the 

18 unknown coefficient vector y and the matrix T . 



-47- 



1 32. The method of claim 31 wherein the matrix 

2 T is based on three sets of samples derived from the 

3 desired pulse shape, wherein the samples have a period T, 

4 wherein the first set comprises the pulse shape samples 

5 shifted by +T/2, wherein the second set comprises the 

6 pulse shape samples not shifted, and wherein the third 

7 set comprises the pulse shape samples shifted by -172. 

1 33 . The method of claim 32 wherein the desired 

2 pulse shape comprises a complex raised cosine pulse 

3 shape . 

1 34. The method of claim 31 wherein the matrix 

2 r is constructed with a number of zeros dependent on 

3 distances between peaks of the correlation. 

1 35. The method of claim 31 wherein the 

2 determining of the channel response from the unknown 

3 coefficient vector y and the matrix T comprises 

4 multiplying the matrix V by the unknown coefficient 

5 vector y . 
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1 36. The method of claim 31 wherein the desired 

2 pulse shape comprises a complex raised cosine pulse 

3 shape . 

1 37. The method of claim 31 wherein the solving 

2 of the least-squares solution comprises solving the 

3 least-squares solution using a singular value 

4 decomposition method. 

1 38. The method of claim 31 wherein the solving 

2 of the least-squares solution comprises solving the 

3 least-squares solution using a conjugate gradient method. 

1 39. The method of claim 31 further comprising 

2 adjusting the tap weights of the taps of the equalizer 

3 based on the channel response. 
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1 40. The method of claim 31 further comprising 

2 determining the initial tap weights for the equalizer 

3 based on the channel response and a noise variance <7? 

4 given by the following equation: 
5 

6 <7?= 1 ||y-y|| 2 

7 

8 wherein D denotes distance between peaks of the 

9 correlation, wherein the superscript a denotes a distance 

10 of a corresponding anti-causal peak, wherein the 

11 superscript c denotes a distance of a corresponding 

12 causal peak, wherein K denotes the number of anti-causal 

13 peaks, wherein L denotes the number of causal peaks, 

14 wherein N denotes a number of training symbols in the 

15 training sequence, and wherein y is based on the vector 

16 y. 
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(57) Abstract: A response of a channel may be estimated by correlating a received signal and a training sequence, by forming a 
matrix T based on a desired shape for the peaks of the correlation, by extracting a vector y from the received signal, and by estimating 
^ the channel response from a least-squares solution based on the matrix I\ the vector y, and a matrix formed from the elements of the 
known training sequence. 
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